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We study the influence of the fluctuations of a Lorentz invariant and conserved vacuum on cos¬ 
mological metric perturbations, and show that they generically blow up in the IR. We compute this 
effect using the Kallen-Lehmann spectral representation of stress correlators in generic quantum 
field theories, as well as the holographic bound on their entanglement entropy, both leading to an 
IR cut-off that scales as the fifth power of the highest UV scale (in Planck units). One may view 
this as analogous to the Heisenberg uncertainty principle, which is imposed on the phase space 
of gravitational theories by the Einstein constraint equations. The leading effect on cosmological 
observables come from anisotropic vacuum stresses which imply: i) any extension of the standard 
model of particle physics can only have masses (or resonances) < 24 TeV, and li) perturbative 
quantum field theory or quantum gravity become strongly coupled beyond a UV scale of A < 1 
PeV. Such a low strong coupling scale is independently motivated by the Higgs hierarchy problem. 
This result, which we dub the cosmological non-constant problem, can be viewed as an extension of 
the cosmological constant (CC) problem, demonstrating the non-trivial UV-IR coupling and (yet 
another) limitation of effective field theory in gravity. However, it is more severe than the old CC 
problem, as vacuum fluctuations cannot be tuned to cancel due to the positivity of spectral densities 
or entropy. We thus predict that future advances in cosmological observations and collider technol¬ 
ogy will sandwich from above and below, and eventually discover, new (non-perturbative) physics 
beyond the Standard Model within the TeV-PeV energy range. 


I. INTRODUCTION 

The cosmological constant problem is arguably one of the deepest and most long-standing puzzles in 
theoretical physics [1, 2]. Naively one would expect energy in the vacuum to couple to the spacetime metric 
and drive accelerated expansion at an enormous rate, governed by a UV scale A above which there are is 
no coupling to gravity from modes in the vacuum. The observation of a large and very slowly accelerating 
universe [3, 4] raises the question of why the vacuum does not gravitate in this way (or why A ~ meV 
is so small, while quantum field theory is well tested up to energies 15 orders of magnitude higher), and 
much effort has been invested in understanding this apparent gross inconsistency [5-7]. At this level, 
the cosmological constant problem asks how the expectation value of the vacuum stress-energy tensor 
affects the metric as a source in Einstein’s equation. Quantum mechanically, the vacuum is 
characterized by higher order correlators, beginning with the two-point function {Tll^\x)T^]^\y)). Just as 
the vacuum expectation value ) can influence the homogeneous background cosmology, these higher- 
order correlations influence higher-order classical statistics of (inhomogeneous) cosmological perturbations. 
Thus, there is a Cosmological non-Constant (or CnC) problem which asks how the ultraviolet physics of 
the vacuum encoded in these correlators affects spacetime geometry on infrared, cosmological scales. 

In particular, we will here address the question of how two-point functions for metric perturbations can 
be affected by the vacuum, and whether cosmological observations can place constraints on UV physics. 
While we find that energy and momentum density fluctuations are small in the infrared, the vacuum stress 
components scale as 


(I;7)(k)^i^\k'))~J3(k + k')A^ 
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leading to point-like correlation in real space, 


( 2 ) 


The gravitational potential $ is affected by anisotropic stress as ~ where is a tensor 

which picks out the anisotropic part of and Mp = is the reduced Planck mass. So Eq. (1) 

yields a contribution 


(A^ 


A5 

Jpk 


( 3 ) 


to the dimensionless power spectrum A|,, defined by 


($k$k') = (27r)353(k + k')A|(fc) (fcVSTT^) ' ■ 


( 4 ) 


These fluctuations will affect observable quantities. To a first-order approximation, we can require that 
A|> ^ 1 to maintain consistency with the observation of a homogeneous background cosmology. Requiring 
that A^'^^ < 1 on the largest accessible scales, k ~ Hq, leads to the condition 

A < ~ 2 PeV. (5) 

In the following sections we will make this schematic statement more precise, and strengthen the bound 
by studying how the UV scale can be constrained with cosmological information about the amplitude of 

$. 

Vacuum stress-energy fluctuations - in particular two-point correlators - have been studied in a fiat 
space context [8] as well as on curved space [9, 10], including de Sitter space [11, 12] and inflation [13]. 
The physical effects of stress tensor fluctuations on the gravitational field were reviewed in [14] . We also 
note that our results - in particular Section VI - share some overlap with results in stochastic gravity (see 
[15] and references therein). 

Stress-energy correlators generically feature UV divergences, requiring a method of regularization or 
renormalization to obtain a finite and physically meaningful result. We will encounter this obstacle in our 
treatment below, and note that regularization of stress-energy correlators has been studied extensively in 
the literature (e.g., see [16, 17] and references therein). 

In Section II, we present Einstein’s constraint equations for scalar and vector perturbations in a fiat 
FRW background, sourced by stress-energy perturbations. We do not consider tensor perturbations in this 
work, as they do not appear in linearized Einstein constraint equations, and thus are much less sensitive 
to the short-wavelength vacuum fluctuations, due to time averaging over the past light cone. For further 
discussion of this point see Section VIII. 

In Section III, we parametrize vacuum stress-energy in a general quantum field theory context using 
the Kallen-Lehmann spectral representation, and identify the large-scale behavior of metric perturbations. 
In Section IV, we study the influence of vacuum stress-energy on gravitational forces on astrophysical 
objects and the Hubble law. In Section V, we introduce a classical toy model for a conserved, generally 
covariant two-point function for stress-energy, featuring a lack of spatial correlation, or Poisson statistics. 
We compute the influence of this matter source on cosmological perturbations and on the cosmic microwave 
background (CMB) through the integrated Sachs-Wolfe (ISW) effect. In Section VI, we study vacuum 
stress-energy for a massive scalar field and find a spectral density closely connected to that of the Poisson 
model of Section V, through analytic continuation. Finally, in Section VII we present an independent 
constraint on UV physics from requiring that the entropy of a region of space not exceed the holographic 
bound given by the Bekenstein-Hawking area law. We discuss our results and conclude in Section VIII. 

Those interested in our results may skip to the ends of sections IV, V, and VII, and note the boxed 
equations there. We will use the (—1"++) metric signature, and natural units = c = fcs = 1 throughout. 
We will denote the scalar product of four-momentum as = fcl to distinguish it from the magnitude 
of three-momentum |k| = k. 
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II. METRIC PERTURBATIONS IN FRW BACKGROUND 


We work in the conformal Newtonian gauge for metric perturbations, 

ds^ = a‘^{r]) [—(1 + 2(j))drf‘ + 2Vidxidri + (1 — 2?/;)dx^] . (6) 

The Einstein equations for scalar metric perturbations are, in Fourier space, 

+ = 4TrGa‘^ST°, (7) 

+ il)') = 47rGa^(5r°, (8) 

(V'" + n{2il; + (j))' + {2n' + + \k,k,{(t) -^) = AnGa^ST;, (9) 


where 0 components and ' derivatives refer to conformal time rj. Eqs. (7) - (9) can be solved for tjj and (j) 
in terms of the stress-energy, 

ik^5T,^ , (10) 

—k'^cj) = A-kG ^i5Too — 




sn 


(11) 


— k^ij} = 47rG ( STqo — 




as well as the time derivatives 


— k'^tp' 
—k‘^(j)' 


AnG 

AtvG 


-H5Tqo + 1 + 


- H6Too + 1 + 


fc2 

377^ 


ifcMT.o - -H - 3 


k^k^ 


fc2 


fc2 

k'^y 




tyST^o + ( 6^^ - 3— ) {6Tk - nST^j) 


( 12 ) 

(13) 


Note that in the subhorizon regime k ^ H, at leading order we can drop 0{'H) and OiTi,^) terms. 


-fcV = IttG^Too, (14) 

-k^cj) = 47rG[(5roo + (<5*^-3fc*F)(5Ty], (15) 

—yip' = AnGiySTiQ^ (16) 

-k^(p' = ATrG[iy6T,o + {S'^-3yy){STijy], (17) 


where ki = ki/k. 

For vector metric perturbations, Einstein’s equations take the form 

yVi = lQ'KG{5ij — kikj)5TjQ, (18) 

where the projection {5ij — kikj) picks out the purely vector part of STjq. The conservation of stress-energy, 
^rpfiu _ gjygg yg Fourier space constraints 

((5r°)' = ik'sry -3nST° + nST^, ( 19 ) 

{ST°y = -Amry - iy 5Ti, ( 20 ) 

from which we see that 

yv' = IGnGiSij - hkj){iy6Tji - 2nSTjo). (21) 

Eqs. (10)-(13), (18), (21) show the impact of vacuum stress-energy on metric perturbations. In subsequent 
sections we will consider specific sources 6T^^. 

At the linear level, the Einstein equations hold unambiguously as equations for Heisenberg operators, 
with metric fluctuations treated quantum mechanically in the same way as stress-energy. 

The effect of vacuum stress-energy on metric fluctuations which we compute is exact for linear per¬ 
turbations around a Minkowski spacetime background. Additional matter, or any form of background 
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stress-energy, will lead to curvature in the geometry. In this case, the effect of vacuum stress-energy on the 
metric backreacts in turn on other (non-vacuum) forms of stress-energy. A fully self-consistent analysis 
would have to account for this in solving Einstein’s equations. The flat-space effect, however, is the leading 
contribution in an expanding universe, with subleading contributions suppressed by "H/Zc on sub-Hubble 
scales. In the super-Hubble regime k <T-L, these corrections become large, and the Minkowski space result 
is no longer the approximate solution. 


III. VACUUM STRESS-ENERGY IN THE KALLEN-LEHMANN REPRESENTATION 


We now move on to consider stress-energy fluctuations in a general quantum field theory context. We 
can parametrize the non- time-ordered two point function for a conserved, Lorentz invariant stress-energy 
source in the Kallen-Lehmann spectral representation [18] as 


c^s — g {To,p(y)) 


drk 


Jk-{x-y) 


dfi 


PGif^)Pfii'Pal3 “f P2 (m) ( 


^Pyti/Pap 


2'K5{k\ 


■p), 

( 22 ) 


where the projection tensors = 77 ^,^ — k^k^/kl ensure stress-energy conservation = 0. Here, the 

c subscript denotes the connected part of the two-point function; in what follows we will omit this subscript 
and implicitly discuss only the connected part. ()s denotes the symmetrized two-point correlation, which 
makes it the expectation value of a Hermitian operator, and thus a potential quantum observable. In 
Appendix A, we discuss the tensor structure and spectral densities po,2, and show how scalar contractions 
of Eq. (22) relate to the spectral densities. The subscripts refer to spin-0 or spin-2 states; as shown in 
Appendix A, only the tensor structure for the spin-0 part contributes to correlators with the stress-energy 
trace T^; the spin-2 tensor structure comes from the traceless part of Eq. (22) is consistent with 
unitarity if and only if / 9 o, 2 (m) > 0 .^ 

It is important to note that for UV physics, only high frequencies will contribute to {T^^{x)Tap{y))', the 
spectral densities po ,2 vanish for y smaller than the masses of heavy particles. 

Fourier transforming the spatial arguments and defining T^,^(k) = J we have at equal 

times 


Jo 






P L 


P2 (/^) ( P 2^PfipPiya. ^PfiipPap 


where in and the equal-time power spectra are defined by 


(T^y(k,t)r^y(k',t)) = (27r)353(k + k')P^2/3(k,i)- 




(23) 


(24) 


Since we are interested in the leading effect on IR scales k = |k| from physics at UV scales fi, we neglect 
k in these expressions; that is, we drop terms at higher orders in k'^ /y. We define 


1 Pf-^max 

= 2 s = 0,2. (25) 

We will drop the argument for fi^nax 00. Note that Jm has mass dimension 2 to -|- 6 , and in particular 
J-ij 2 , which is most relevant here, is dimension five. 

Due to the projection tensors in Eq. (22), which impose stress-energy conservation, long-wavelength 
vacuum pressure and stress fluctuations are large on infrared scales, and energy and momentum density 
fluctuations are suppressed in comparison. This is because 0 indices lead to a suppression by k/ so 


^ This may be shown by requiring that the expectation value {T{J)^) of the square of T{J) = f be non-negative 

for an arbitrary real-valued source function Jk-''(x). We thank Andrew Tolley for correspondence on this point. 
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that only spatial components (TijTki) are dominant. At leading order, then, 

t) = + (^ls.k6,i + J%^ + Oik^), (26) 

and components of (T^j^^Tap) with fewer spatial indices are suppressed by powers of k. 


A. Cosmological Influence of Vacuum Stress-Energy 


We now treat the stress-energy of Eq. (22) as a perturbation in Einstein’s equations, and find 

the contribution to metric (equal-time) two-point functions sourced by the stress-energy as given in Eq. 
(22). We will place a V superscript on power spectra to indicate that only the contribution to two-point 
functions sourced by the vacuum are shown, not the full two-point functions. The “0” index will indicate 
conformal time r). 

In Einstein’s equations, (11)-(13), the traceless tensor (d®-' — iPk^) picks out only anisotropic stress 
to source metric perturbations, and the transverse projection for vector modes in Eq. (21) has the same 
effect. Thus, (isotropic) pressure from the vacuum does not couple to metric perturbations, and as noted 
above, other components of are small in the infrared, leaving only anisotropic stress as the vacuum 
source of infrared geometry. We consider the subhorizon regime fc 3> H. Going from a Minkowski to FRW 
geometry^ introduces a scale factor, since k^P^{k) is a physical quantity (it sets the real-space variance of 
(j)) and should therefore depend only on the physical wavenumber, ~ 1/fcph = a/k. The same applies to 
k^P<p'(k)/, with an additional coming from converting physical time d/dt to comoving time d/dij. 

Using Eq. (26) with Eqs. (15) and (17), then, we find 


k^P'^PiKa) 

lfc3p(r)(fc,a) 


3 '^- 1/2 a 

2 k 

3 '®i/2 a _ 
2l^k ~ 


3a 

" 4M4fc 
3a 

4M4fc j 


I dlJLy/Ji,p2{lJi) 


(27) 

(28) 


with P^'^ and P^ij{k) suppressed by 0{k/,/Ji). In the second line, the time derivative acting on Eq. (22) 
simply adds a factor of fcg « —= /r in the integral. Here, Mp = (SttG)”^/^ is the reduced Planck mass. 

In the following section, we will use this result to constrain the spectral density p 2 via its effect on 
gravitational forces ~ V(j) on astrophysical objects. Although only high frequencies contribute to the 
integrals in Eqs. (27)-(28), we argue in Section VI below that for a finite range of time, transforming 
stress-energy correlators to real time will lead to correlations at large time separation. Furthermore, two 
high frequency modes contributing to (j) can generate at nonlinear order a low frequency mode. We consider 
a resulting cosmological effect in the following section. 


IV. CONSTRAINTS FROM THE HUBBLE FLOW 


Here we consider the gravitational force V0(x) on galaxy clusters in the kinematic Sunyaev-Zel’dovich 
effect, and identify constraints on vacuum stress-energy via its contribution to this force. 

Take the action for a relativistic point particle. 


Sp = —771 


2(j){x,t) - 


(29) 


where, following the discussion above, we neglect other metric perturbations compared to 0, so that 


^ We also expect the curvature to introduce contributions scaling as "H/Zc, which would likely prove difficult to compute, 
since we would be computing vacuum expectation values in curved space, where there is no unique vacuum. In Section 
VIII we comment on the case of (Anti) de Sitter space, where we may be able to compute vacuum stress-energy correlators. 
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~goo = 1 + 20. Furthermore, the proper time of an observer at the origin is given by: 


T = 


+ 20(O,t). 


( 30 ) 


Now, expressing the particle action to second order in 0 and x, in terms of the observer’s proper time, 
yields 


50 ~ m 


dr 


-1 + + 0 ( 0 , t) - 0 (x,t) + ^ 0 (x,t )2 - ^ 0 (O,t)^ + 0 (x,t) 0 (O,t) 


(31) 


Given that 0 is oscillating at high frequencies (given by the UV scale of the theory), the linear terms cannot 
affect the dynamics of macroscopic observables (e.g. a galaxy). In fact, averaging over the oscillations, the 
only non-vanishing and non-constant contribution to the effective Newtonian potential is the cross term: 

$w(x,f) ~ -(0(x,f)0(O,t))high, (32) 

where we have taken the quantum expectation value only over high frequencies. Averaging over low 
frequencies gives 


/ d?k 

exp(jk • x)P 0 (fc), (33) 

where P^{k) oc k~^ is given by Eq. (27). While the integral in (or the 0-correlation function) is 
formally divergent, the acceleration is finite and given by: 

( 2 ) 

(x-V$Ar(x,t))--y yy^(ik-x)exp(ik-x)P,^(fc) = ^^^|x|. (34) 

The surprising conclusion is that, at second order in perturbations, the high frequency oscillations lead 
to a uniform and constant centripetal acceleration. While Eq. (34) provides the perturbation to the 
Newtonian gravitational force due to high frequency fluctuations, its integral over the classical trajectory 
provides the change in radial velocity. Within linear perturbation theory, this leads to a constant infall 
peculiar velocity, or a negative offset to the Hubble law:^ 


V = i7x — 


gi^m) 

H 


V 


= iJx- 


■^-1/2 X 

167 riJ |x| ’ 


(35) 


where comes from solving for radial trajectories, to linear order in perturbation theory in a ACDM 

background: 


giS^rn) — 


/p (1 - a^)h{a,Vtmy 


/g h{a, Q.m) ^da 


I ! - 

-, flm) = V ^ -f 1 fin 


(36) 


By measuring the kinematic Sunyaev-Zel’dovich effect of X-ray selected galaxy clusters, the Planck 
satellite [19] has measured the mean radial peculiar velocity to be {vr) = 72 ± 60 km/s, which implies 
(vr) > —25 km/s = —8.3 x 10“®c at the 95% confidence level. Plugging this into Eq. (35), and using the 
( 2 ) 

definition of (^ 9 - yields: 


t(2) 

[^-1/2 

1/5 


y/y \ 


-I 1/5 


< 1.0 PeV, 


( 37 ) 


^ Observers elsewhere would see the same effect in their proper time, which depends differently on 4 >{'x.). 
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at the 95% confidence level. 

Given that / 02 (m) at high energies for a weakly coupled quantum field theory (e.g., see Eq. 72 

below), we can interpret this upper limit as an upper limit on the energy scale of perturbative quantum 
field theory, or quantum gravity. Beyond this scale, the theory must be strongly coupled (with no weakly 
coupled local UV completion), as /92 (m) must decay faster than for the J^i /2 integral to converge. 


V. A POISSON SPRINKLING OF PHASE SPACE 

In this section, we consider a toy model for the vacuum comprised of a classical collection of massive 
particles in phase space, which is manifestly Lorentz invariant and conserved. We will see that a Lorentz 
invariant stress-energy correlator describes this model, along with a fiat Poisson spectrum, as in Eq. (1), 
and will find an ISW plateau for masses > 24 TeV, which is already constrained by CMB observations. 
We shall show in Section VI below that a free massive scalar quantum field theory reproduces an identical 
structure to this toy model on large scales. Here and in the following subsection we restrict to Minkowski 
space; in Section VB we generalize to an FRW background. 

With a phase-space distribution /, the stress-energy tensor integrates over momenta with the Lorentz 
invariant measure 


T^''(x,t) = y ^p>'^/(x,p,t), (38) 

Here, is the four-momentum for classical point particles of mass m. We consider a distribution /(x, p, t) 
characterized (like Eq. (2)) by point-like correlation or Poisson noise^ in phase space, 

(/(x,p,<)/(x',p',t)) - (/(x,p,t))(/(x',p',t)) = {f{x,p,t))S^{x-x')S^{p-p'). (39) 

The expectation value here is taken with respect to an ensemble of distributions of particles, from which 
a particular realization / is drawn. From now on we will ignore the constant (/)^ term, as this will only 
shift the background, and correlators for stress-energy fluctuations. We take the mean distribution to be 
uniform in phase space®, 


(/(x,p,t)) = (/(t)). (40) 

We will see below that this results in a Lorentz invariant two-point function. 

To evaluate the stress-energy two-point function, we note first that the phase space distribution is related 
to the initial distribution /o at time to by 


/(x.p,t) = /o 



{t - to) 


dx 

dt 



(41) 


Putting together Eqs. (39) and (41), we find 


(r^‘'(x, + Ax,t + At))c 


d^p 

w 


P»V(/o)5" {Ax^ 


AtpV/) . 


(42) 


For spacelike separation |Ax| > At, the delta function vanishes for all p. For 


timelike separation it fixes 


Ax^ 

p^ = m — , —, 

-ly—AxyAx'^ 


(43) 


and contributes a normalization factor m^At^{—AxjAx'^) That is, the only particles in the original 
collection contributing to the correlation function at separation Ax^ are those with the corresponding 


For zero correlation between different phase space points as in Eq. (39), the probability for finding N particles in a given 
region of phase space is a Poisson distribution. 

^ To keep the distribution normalizable, we assume that at high momenta p ^ m, (/(p)) falls to zero. 






four-momentum given by Eq. (43). Putting everything together, we have 


{T^''{x)T-P{x + Ax)), = e(-Ax^Ax^). (44) 

Eq. (44) is a conserved, generally covariant two-point function in Minkowski space. We will use it as a 
proxy for vacuum stress-energy, and argue below that the two behave similarly on large scales. 


A. Spectral Density for Poisson Stress-Energy 


The spectral densities for Poisson stress-energy in the Kallen-Lehmann representation cair be found by 
transforming to Fourier space. Defining T^,j{k) = f fi'^xe“*^'^T^,y(x), we have 


{T^u{k)T^p{k')) = {2n)HHk + k') J + y)). 


(45) 


Using Eq. (A5) to relate the left hand side to po{—kl), we find 

27Tp^~(-k!) = J d^xe^^-^T^,{y)T^p{x + y)),.,. (46) 

For Poisson stress-energy the Fourier transform can be most easily obtained from Eq. (42). Eliminating 
the spatial integral with the delta function there, and integrating over time, we find 

/Oo°i""°n(-fc|) = ^{fo)Tn‘^ J ^6{k^pf^) = J d'^pSip'^ + m'^)6{kf,p>")Q{po). (47) 

For fcl < 0, vanishes due to the delta function since p^ < 0. For k\ > 0, we can evaluate the 

integral in a Lorentz invariant manner by replacing the delta functions with integrals over phases, and 
using the stationary phase approximation: 

[ d^p6{p^ +m^)5ik^p>^)Q{po)= r (48) 

J J — oo d pq>0 

The phase is stationary with respect to p at 2xp^ -I- = 0, so 



where in the first line the 7r^/2x^ comes from evaluating the p integrals after shifting p —>■ —ky/2x 
and extracting the phase, and in the secoird line we have integrated over y. The integral J dze''^ 
is divergent. However, we can deform the coirtour to integrate along down the positive imaginary axis, 
around the origin, aird back up with a small positive real part, so that J dze^^ j —>■ ^ dze~^ jz^!"^ = 

r(-i/2)M = —2^Tili. The appearance of factors of y/i reveal an ambiguity in our regularization scheme; 
we choose the prescription which gives positive spectral densities (see footnote 1). 

Also recall that since the Poisson stress-energy correlator is totally symmetric in its indices, its tensor 
structure in the Kallen-Lehmann representation will be totally symmetric, Pfi^Pap + Pp.aPvp + Pp.pPiya, 
which fixes 


Poisson //z /c:^ Poisson 

P2 = (6/5)Po 


( 50 ) 
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Putting everything together, then, we have 


Poisson 

F 2 


i-kl) 



(51) 


In the following section we will make use of this result to study the cosmological influence of Poisson 
stress-energy. 

The regularization of the divergence at small x or z in Eq. (49) controls the divergent behavior of Poisson 
stress-energy at large four-momentum p, which corresponds to approaching the light cone Ax^ = 0 , as seen 
in Eq. (43). This gives us a well-behaved result for power spectra which, matches the IR limit of quantum 
field theory stress-energy correlators in the Kallen-Lehmann representation, as we will see in section VI. 


B. Cosmological Influence of Poisson Stress-Energy 


We can now And the contribution to metric two-point functions sourced by the stress-energy as given in 
Eq. (44), which we take to be a small perturbation . We will place a (P) superscript on power spectra 
to indicate that only the contribution to two-point functions sourced by Poisson stress-energy is shown, 
not the full two-point functions. Making use of the symmetrical index structure of Poisson stress-energy, 
which leads to Eq. (50), we have 

{6T^lKk)ST^^^Hk')), = {27:fSHk + k')27Tp^~{-kl)^ + P>.o.P.p + • (52) 

For the Poisson spectral density, Eq. (51), we And from Einstein’s equations in the fc = |k| ^ H regime 
the following results for two-point correlators in momentum space, where we define (V(k, w)V(k',a;')) = 
( 27 r)^( 5 ^(k -I- k')^(a; -I- uj')Px{k,uj), 



PPik,u;) 

P^%ik,u;) 


9 ^ 

TT am 
10~M^ 


ifo) 


CJ 


- a ;2 


1-1 


0 (fc^ — w^). 


40 M4 
27r^ am^ 

Is" 


ifo) 


(fc2 _ ^2)5/2 


0 (fc^ — 


{fo){Sij - kikj) 


—OJ 


fc2 {k^-OJ^fA 


0 (fc^ — W^). 


(53) 

(54) 

(55) 


We have included the scale factor a{r]) to account for the effect of FRW space, through the follow¬ 
ing argument: While the power spectrum P^i{k,ui) is not a physical quantity itself, the combination 
k^ujP^i{k,uj)/a^ is a physical quantity, since the real-space correlation ((/)(xi,<i)(()(x 2 ,^ 2 )) for a small 
range of momenta (Ak, Aw) peaked around (k, w) goes like this combination. Thus, k^ujPtp'{k,uj)/a? 
should depend only on physical frequencies and wavenumbers, fcph = k/a and Wph = w/a. Making this 
replacement in Eqs. (53)-(55) adds the factor of 0 ( 77 ), as shown above. A time-dependence to the frequency 
spectrum is valid for sub-Hubble times oj > H ( 77 ). 

Upon integrating over frequencies w to obtain correlators in real time, we obtain divergent pieces going 
as powers of Wmax, but for P^,' and Py, a finite, positive part remains (despite the negativity of Py,^), 
which we take to be the physical effect of vacuum stress-energy. We also note that as emphasized in section 
III A this effect comes from the correlator of spatial components (STijSTki), and from anisotropic stress, 
that is from p 2 rather than the isotropic po part of Eq. (22). 


C. ISW Effect from Poisson Stress-Energy. 

We now show that Poisson stress-energy contributes a flat spectrum (independent of 1) to the CMB 
angular power spectrum, contributing most significantly at high I through the ISW effect. Including the 
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influence of vector modes [20] , the ISW effect is 

/*^today 

— = / dil{4'' + V'' + ylr"), (56) 

•'Vlss 

Here, the integration is carried out over the past light cone, that is (j)' = 4>'{r,r,r]) with r{r])/rLss « 
1 — ’q/rjtoAay The contribution to the CMB angular power spectrum is well approximated by the Limber 
approximation [21], which, for integration over the light cone, gives 


5risw( 


qISW 


c?H dr 
Att 





(P) 


V',ij 



(57) 


where fcj_ = (/ + \)/r. Here, we have dropped the (V’VO power spectrum and cross-spectrum {cj/ij}'), 
which only have power-law divergent contributions after integration over w, as well as the cross-spectra 
{(j)'Vl) and (^'17'), which vanish due to the transverse nature of the vector modes. We have also dropped 
power-law divergent contributions to the {4i'4>') {VlV^) correlators. Furthermore, for the vector term, 
only the isotropic part of Eq. (55) contributes. 

The power spectra in Eqs. (53) and (55) go like l/fc^^ after integration over w, canceling the r-dependence 
in the integrand of Eq. (57). However, the scale factor in Eqs. (53)-(55) leads to the radial integral 
contributing / dra{ri^ss ~ ’’) ~ / dr]a{r]) = J dt ^ for = 0.3 and Ha = 0.7. Making use of Eqs. 
(53) and (55) and using (H^jHa) = (0.3, 0.7), we And that at high I 


. a 2 nISW ^ ?(^ + l)Cf^ ^ 497r mHo ^ 497r 

^ ~ 2 tt 720 “ 720 M^iJo ’ 


(58) 


which is a flat spectrum with Z, and thus largest compared to CMB anisotropies in linear theory at the 
end of the damping tail at high 1. In Appendix B, we show that at the 95% confidence level, current CMB 
observations roughly translate to: 


Since Mp = 2.44 x 10^® GeV and Hq « 70 km/s/Mpc = 1.5 x 10“®® eV, requiring the ISW contribution, 
Eq. (58), to be consistent with the observed values implies 


m < (24 TeV) 


r(2^)"(/o)i 

-1/5 

2 



(60) 


In Section VI, we motivate this choice for (/o) by drawing an analogy to quantum held theory of a free 
massive scalar held. 


VI. POISSON-LIKE STRESS-ENERGY FOR A MASSIVE SCALAR FIELD 

In this section, we compute the stress-energy of a massive scalar field, and show that it behaves similarly 
to the Poisson spectrum considered above. While we And that due to the requirement of positive energy, 
the amplitude of stress-energy fluctuations vanishes for small |fc^|, for large enough time separations At 
the spectral density at high energies contributes to the stress-energy correlator in the same way as the 
low-energy Poisson spectral density. 

For a massive scalar field we have 


[(^v?)^ + . (61) 

Defining the momentum space quantities 

J d‘^xe~^'''^(p{x), T^^{k) = j d'^xe~^'^''^Ti^^{x), 


= 


( 62 ) 















we have 
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= j k - P)PpPk-p, (63) 

where 

k') = -k^kl + ]^gtiiy{k°‘k'^ - rr?). (64) 

The two-point correlation is then given by 

/ ^4 ^4 / 

- p)M^/3{p', k' - p'){ipp<fk-pPp"Pk'-p')- (65) 

The scalar field two-point correlation can be evaluated in terms of commutation relations for creation and 
annihilation operators: 

{(fikPk') = {2Tr)‘^6'^{k + k')2'K5{k^ + TO^)0(fco). (66) 

Evaluating the connected four-point correlation in the same way, we find 

{^p^ipk-pPp'Pk'-p')c = (2T^)^8^{k + k')8{p^+rn^)[2'Kf8{{k--p)^^vr?) 

X [(27r)^(5'^(p -f p') -h (27r)'^5^(A: - p -f fc')]0(po)0(A:o - Po)- (67) 


Consequently, letting p —>■ —p in the integral, we have 

{T^„{k)Tai3{k'))c = S'^ik + k') J d'^pM^^{p,k - p)Map{p,k - p) 

X (27r)^(5(p^ -I- m‘^)6{{k - p)^ -I- m^)0(po)0(fco - Po)- 
To find p 2 via Eq. (A6), we contract with , 

S' 2 ‘"“^M^i,(p, k - p)Map{p, ^ “ P) = ^ + 2m^fc^ -I- 8m^ -I- (2fc^ -I- 4m^)p • fc -I- 2(p • fc)^] . 


( 68 ) 


Here we have used the fact that p^ = —We also have 


k - p)Map{p,k - p) = ^[k'^ -p-kY . 

For timelike k we are free to choose k = 0, so fc • p = —kopQ, and 

2 p 2 (-fc^) 0 (fco) = / 7^^(27r)(5(p^-f TO^)5((fc-p)2-hm^)0(po)0 (fco-Po) 


(69) 


(70) 




(27r) 

I ^ sJklH-mPB[k„ - 2m)^ 


21 2 

1 - 4 — 

h'^ ' 

'^0 J 


(71) 


where Wp = -y/lpp -I- m?. In the second line, the delta function sets the step function equal to unity, but 
requires the addition of the step function in the third line. The factor of 20(fco) indicates the unsym¬ 
metrized expectation value. Symmetrizing as in Eq. (A6), and replacing k^ —>■ —fc| to write the expression 
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in its covariant form, we have 


P2{-kl) = 


1 


fcl 


240 167r2 


to 

r 21 


^4 


0(—^4 — 4 m ^). 


(72) 


(For spacelike k on the other hand, we are free to set ko = 0, leading to Q{po)Q{—po) = 0 and a vanishing 
correlator, so that Eq. (72) applies for all k.) In the same way, we find from Eq. (70) that 


Po{-kl) 


1 ki 

18 167r2 



1 

2 



e{-kl - 


4m^). 


(73) 


Eqs. (72)-(73) are consistent with existing literature on scalar-field stress-energy; see e.g. [15, 22], 

In order to find stress-energy correlators in real time at a separation At, we integrate over frequencies 
kg, as in Eqs. (22)-(23). Physically, in the absence of very energetic particles we do not expect linear 
macroscopic observables to be sensitive to very high frequency fluctuations in the metric. If we lack access 
to energies above some scale uJma.x, the stress-energy correlators that are physically relevant for metric 
perturbations will be approximately given by a Fourier transform over frequencies that is damped above 

^max 


/ oo 

dkoP2(k^ - . (74) 

-OO 

For |A;|| ^ becomes a gaussian integral, which can be evaluated analytically, and is equal to a poly¬ 

nomial X exp (—This allows for the integration over kg, which runs over |fco| > -s/k^ + 4m2, 
to be done along the contour Coo shown in Figure 1, as long as it is closed at \k1\ ^ rn?. The parts of 
the contour around the branch cuts along the real axis give us four times the original integral, while the 
integration at |fco| —t oo can be dropped for cjinax|At| 1. We can then shrink the contour to Ci, which 
encloses the branch points at kg = ±|k|. Taking Wmax |k| and dropping terms suppressed by [kp/m^ 
in Eq. (72), this is 


Jp{k, At) = 


m 


1207r2 


dkg 


-fcg -I- k 2 


ifcoAi 


for |At| uj. 


-1 

max' 


(75) 


So the scalar-field stress-energy contributes in the same way as a spectral density at low frequencies 


effective 

r2,ip 


i-ki) 


m 


1207r2 



(76) 


integrated for real kg. 

The scalar field spectral density, Eq. (72) , therefore contributes to correlations in real time in the same 
way as the effective spectral density, Eq. (76). This is equivalent to the low-frequency spectral density for 
the Poisson stress-energy of Section V A, so we expect the same stress-energy correlations at long times At 
and large distances |k| <C m which led to the TeV scale bound discussed there. (The Poisson stress-energy 
correlation is that of a collection of pointlike particles, which we only ought to compare to the quantum 
field theory case on spatial scales much larger than the Compton wavelength m~^ of the particles.) 

Comparing to Eq. (76), we see that 


\P 


effective 

2,if 


_ _ _ Poisson 


(77) 


Motivated by the similar behavior of these spectral densities, we consider equating them to fix (fg) = 
2/(27r)^. As noted in Section VC, this places the bound at m < 24 TeV. We note that the ratio of 
the effective spectral densities pg and p2 for a scalar field is equal to the ratio for Poisson stress-energy, 
strengthening the analogy. Following the above procedure for Eq. (73) to pick out the leading term in the 
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FIG. 1: The integrations over Coo and Ci are equivalent. 


k1 limit leads to 


effective/ / 2\ effective/ i 2\ /'to\ 

(-^ 4 ) = K-ki)- (78) 

We conjecture that the same Poisson behavior arises generically in the fc| —>■ 0 limit for massive quantum 
field theories. 


VII. HOLOGRAPHIC ENTANGLEMENT ENTROPY BOUND AND AN IR CUT-OFF OF 

GRAVITY 


In this section, we argue that an IR cut-off in quantum gravity, similar to what we have discovered in the 
CnC problem above, should have also been anticipated from a holographic bound on entanglement entropy 
of non-gravitational degrees of freedom. This independent line of reasoning reinforces the interpretation 
of the scale where metric fluctuations become large (due to stress-energy sources) as a minimum energy 
scale or maximum distance scale for the theory. It also shows the uniqueness of our results to gravity in 
comparison to other gauge theories such as QED (see section VIII D), which lack a holographic entropy 
bound. 

Imagine a quantum system, described by the Hamiltonian Hq in its ground state |0)o, whose energy is 
set to zero for convenience. Now, if we turn on a perturbation Hint, then the static wavefunction |0)o is 
no longer a solution of the Schrodinger equation, but rather will slowly evolve in the Hilbert space of Hq 
with the amplitudes that are (to first order) given by 


(n|0)o ^ 


(n|Hint|0) 


(79) 


where |n) for n > 0 are the excited states of Hq. Even though, at any time, |0)o is still a pure state, we 
can define a density matrix for observation at random times (or averaged over time): 

Pint = ^ |(n|0)onn)(n| (80) 

n n 

which, e.g. for a two-level system, yields a von Neumann entropy 

Squbit = -tr{pint In Pint) ~ Q: [1 - ln(a)] -I- C>(q!^), 


(81) 
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where 


_ l(l|gint|0)P 
El 


(82) 


can be interpreted as the fine structure constant for the interaction. 

Now, let us apply this result to the ground state of a quantum field theory (QFT), as we turn on gravity. 
Of course, one may wonder whether you could start in the global ground state of QFT+gravity. However, 
it is not clear whether such a state even exists, and even if it does, our evolving universe is certainly not 
in its global ground state. Alternatively, we may consider this entropy to be the entanglement entropy 
of the mixed state of QFT in the (pure) ground state of QFT+gravity, which has a similar parametric 
dependence on a. 

Therefore, starting from the QFT ground state, we expect a minimum (entanglement) entropy of 
—aa ln(Q;G) for every qubit of the QFT, where aa ~ E'^/M^ is the gravitational fine structure con¬ 
stant. The minimum number of qubits in a volume, up to an energy scale A, e.g. for a Dirac field is given 

by 


^ = 2 X 2 X Volume x 



2A3 

Stt^ 


X Volume, 


(83) 


where the factors of 2 account for spins up and down, as well as particles and anti-particles. The condition 
that this entropy should not exceed the holographic (Bekenstein-Hawking) entropy implies: 


Sbh = 27rMp X Area > S' = # x ckg [1 — In(aG)] 
For a spherical region of radius R, this yields: 

R ^ Smax 


2A5 [1 + ln(M2/A2)] 


37r2M2 


X Volume. 




A5 [1 + ln(M2/A2)] ’ 


(84) 


(85) 


which, as in previous sections, provides an IR cut-off (or scale of strong coupling) in terms of the UV 
cut-off of the QFT, when gravity is turned on: 


TT A5 [1 + ln(M2/A2); 

^ iW 3^2Af4 ■ 

Requiring that our observable universe fits within this IR cut-off yields: 

Air < Hq — 

=> 


9.5 X IQ-^^eV 
A < 2.4 PeV. 


( 86 ) 


(87) 


While this is not quite as strong as previous bounds, we should note that it does not use any precision 
cosmology data. In fact, one could have obtained this upper limit, nearly a century ago, upon the discovery 
of quantum mechanics and cosmic expansion! 


VIII. DISCUSSION 

The nonlinearity of general relativity couples modes on very different scales. Classically, cosmological 
metric and matter density perturbations evolve nonlinearly and become coupled on very different scales, 
leading to nonlinear structure formation. Quantum mechanically, metric perturbations on cosmological 
scales can also be influenced by short-scale vacuum modes, allowing UV physics to influence classical 
geometry in the IR. 
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A. Summary 

At large distances, we saw that stress-energy fluctuations are largest from the spatial components 
and behave like Poisson noise, being uncorrelated in real space (or, having a flat spectrum in k-space, 
as in Eq. (26)). Using these fluctuations as a source for metric perturbations, we found constraints on 
vacuum stress-energy from measurements of galaxy cluster peculiar velocities (Eq. (37)), which probe the 
gravitational potential. Requiring that metric perturbations remain perturbative on the largest accessible 
scales leads to a similar constraint from Eq. (27). These constraints place bounds on the integrated stress- 
energy spectral density, f given explicitly in terms of the stress-energy tensor in Appendix 

A. We can interpret this as an upper bound on an energy cut-off beyond which quantum held theory, or 
quantum gravity, is strongly coupled. For the Poisson stress-energy correlator, given in Eq. (44), which 
resembles the stress correlators of a massive scalar quantum held theory in the IR, we found even stronger 
constraints on the UV scale through the ISW contribution to the CMB (Eq. (60)). 

The proximity of our upper bound of 24 TeV-1 PeV to the electroweak scale suggests a connection to a 
possible resolution to the Higgs hierarchy problem, which independently requires new physics close to the 
TeV scale. In particular, the CnC problem points to a gravitational resolution to the hierarchy problem 
such as (supersymmetric) large extra dimensions with low scale quantum gravity (e.g., [23, 24]). The 
other possibility is a transition to strongly coupled non-perturbative dynamics, such as technicolor [25], or 
theories with emergent Lorentz symmetry [26] . We thus predict that future advances in collider technology 
and cosmological observations will sandwich, and eventually discover new TeV scale physics. 


B. Directions for Future Work 

Several extensions of our work are possible. A complete account of the influence of the vacuum would 
include backreaction: matter will move on geodesics now modified by the vacuum-influenced geometry, 
leading to a different matter stress-energy which will in turn affect the geometry. We leave a full self- 
consistent analysis to future work. In particular, the influence of the vacuum on geometry and matter 
geodesics on the scales of gravitationally bound systems would be worth investigating. 

While we have focused on the impact of vacuum stress-energy on infrared cosmological modes, we 
have not considered the effect of modes in the super-Hubble regime fc <C 77. In this regime the leading 
effect of metric perturbations is a shift to the locally observed spatial curvature [27]. Our results here 
suggest that metric fluctuations become large in the IR, ~ 1/k, possibly leading to the breakdown 
of FRW as a background spacetime on very large scales. Given a UV scale which sets the amplitude of 
vacuum stress-energy, Einstein’s equations determine an IR scale at which metric fluctuations sourced by 
vacuum stress-energy become large. It would be interesting to seek a nonperturbative understanding of 
the influence of vacuum stress-energy in the fc —>■ 0 limit, and the consequential effect on the background 
geometry in Hubble-sized volumes. 

In particular, the IR divergences encountered here may be well-behaved in de Sitter or Anti de Sitter 
space, with the curvature radius controlling the IR behavior, and a unique vacuum fixed by the symmetry 
of the spacetime, making curved-space QFT calculations possible. It would be interesting to carry out the 
same exercise of finding vacuum stress-energy correlators, and their effect on fluctuations in the geometry, 
in (A)dS, with an exact computation valid on super-curvature scales. The strong coupling on large scales 
may be dual to a confinement in the boundary QFT. 

We have also refrained from discussing the conversion of quantum fluctuations of the metric (super¬ 
positions of different gravitational held configurations) into smooth, classical perturbations. In order for 
this to take place, a mechanism of decoherence is needed to effectively measure the spacetime geometry. 
We expect that an environment of particles with some characteristic energy will be sensitive to metric 
fluctuations up to frequencies of the same order, and will select out a classical configuration for metric 
perturbations with lower frequencies. Energies in the primordial universe, during a reheating or hot big 
bang era, will reach a maximum before redshifting, and will thus be most sensitive to vacuum stress-energy 
during this epoch. It would be interesting to apply our results in this context, which would require a bet¬ 
ter understanding of stress-energy as a source for metric perturbations on super-Hubble modes, as noted 
above. 
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C 


A 


D 
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FIG. 2: One-loop radiative corrections to the graviton propagator, in QFT+GR. The solid and dotted lines 

denote matter and graviton propagators respectively. 


C. Renormalization 

Here, we argue that the standard renormalization of UV divergences in QFT does not affect the UV-IR 
coupling that we have discussed in this paper. 

The renormalization program in QFT (e.g., see [18]) accounts for the fact that the bare conpling con¬ 
stants that appear in the action are not directly observable. As snch, the dependence on the UV cnt-off 
that appears in radiative (or loop) quantum corrections to scattering amplitndes or propagators can be 
absorbed in a renormalization of the bare conpling constants of the theory. 

The radiative corrections to the graviton propagator, at one-loop order, are shown in Fig. (2). Diagram 
A is the connter-term associated with the renormalization of the bare cosmological constant and Planck 
mass. Diagrams D and E are standard one-loop corrections to general relativity, which e.g., lead to higher 
order corrections to the Newtonian potential (e.g., [28]). Diagram C has the same exact strnctnre as A, 
and thns can be absorbed into the renormalization of the Planck mass and cosmological constant (np to 
natnralness considerations). 

Diagram B is exactly what we have compnted in Section VI, and is the only non-trivial one-loop 
correction to the connected part of the metric two-point correlation function, dne its interaction with a 
massive field. Fnrthermore, we saw that it is finite, and thns does not require renormalization. Another 
way to see this is to note that tree-level counter-terms renormalize the one-loop radiative corrections. 
However, in GR, scalar and vector modes are not dynamical, and thus their propagator vanishes at tree 
level. Therefore, the contribution of Diagram B to scalar and vector two-point correlation functions, which 
we have focnsed on in this paper, cannot be renormalized by GR local counter-terms (or otherwise require 
introducing non-local terms into the GR action). 


D. UV-IR Coupling in Gravity, Electromagnetism, and Chromodynamics 

Our study of the effect of stress-energy sources on the gravitational potential can be compared to the 
analogous case in quantnm electrodynamics (QED), where charge and current density source the vector 
potential in Maxwell’s equation. The conservation law 9^ = 0 again constrains the two-point sonrce 
correlator to have a tensor structnre oc 77 ^^ — in momentnm space, snppressing charge 

density fluctuations in the IR. Current flnctnations may have a flat Poisson spectrnm, (Ji(k) Jj(—k)) c>c Sij, 
which acts as a source to the wave equation for the vector potential Ai in the Gonlomb gange (S^A® = 0). 
This is analogons to what happens to the tensor modes in GR, which also respond to their sonrce through 
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a hyperbolic wave equation. In either case, in contrast to scalar and vector modes that satisfy an elliptic 
equation, the response is an integral over the 4-volume of the past light-cone, rather than a spatial 
3-hypersurface. This leads to additional cancellations due to fast oscillations in time and transverse 
(traceless) conditions, bringing down the degree of divergence from A®/fc to fc^log(A) for tensor modes, 
or log (A) for a vector potential. While the latter is already included in the standard treatments of 
divergences in QED (e.g., [18]), the CnC problem indicates a richer structure of IR divergence in massless 
spin-2 interacting theories, which has been so far overlooked. The infrared divergences of scalar and vector 
metric perturbations studied here are also different from infrared divergences due to soft gravitons in 
Feynman diagrams. The second is a dynamical effect which, like IR divergences in QED, is consistent 
with finite cross sections and observable rates [18, 29[, while the first is a non-dynamical consequence of 
the constraint equations, which cannot be treated in the same way. 

A closer analogy to the CnC problem might be the confinement in quantum chromodynamics (QCD). As 
in the case of the CnC problem, the quantum fluctuations of the gluon field around its free vacuum blow 
up at an IR scale of Aqcd- Below this scale, the gluon field becomes strongly coupled, and an effective 
weakly coupled description is only possible in terms of different degrees of freedom (i.e. hadrons). In the 
case of gravity, this may suggest a non-geometric description on very large scales. However, we caution 
that this analogy must not be pushed too far, as the QCD IR divergence is still logarithmic, and not a 
power law as in the CnC (or CC) problem. 


E. The CnC Problem, CC Problem, and Effective Field Theory 


We emphasize that the cosmological non-constant problem as described here is complementary to, but 
distinct from, the cosmological constant problem. In both cases, a deeper understanding of UV physics is 
needed to explain the absence of any influence on geometry from the vacuum at high energies. However, 
different sign contributions to both contribute positively to the two-point function, so vacuum-sourced 
fluctuations in geometry cannot cancel due to different sign contributions. 

One may wonder why we could not just integrate out heavy fields in the IR, as is commonly done to find 
an effective field theory (EFT). Then, the lore is that UV physics will only renormalize the cosmological 
and Newton’s constants, and otherwise has no observable consequence for low energy observables. In 
contrast, we have clearly demonstrated that the CnC problem has distinct consequences for cosmological 
observations that cannot be mimicked by a change in the cosmological or Newton’s constants. 

To understand this better, we can look at the path integral description of (low energy) QFT+gravity: 

J DgDif X Diff“^[ 5 ,(/j] x exp J d‘^Xy/^{R[g] -f[(/?, g]. (88) 


Here, Diff [g, (/?] is the measure of the path integral that effectively mods out the gauge degrees of freedom, 
in particular diffeomorphisms, in a relativistic theory. What is often called the gravitational EFT is given 
by integrating out the heavy matter fields, but ignoring the measure: 


exp{iSeS , naive b]) = exp(iS'GRb]) X 


Dip exp i / d'^x^/^Cm[^,g] 


(89) 


which can be computed using, e.g., the heat kernel method as an expansion in powers of curvature invariants 
and their gradients (e.g., [30[). However, in general 


Diff ^b,0]exp(iS'eff,„aiveb]) ^ exp(iS'GRb]) X J Dpx Diff '^[g,ip\ x exp j ■ (90) 

Therefore, the standard definition of EFT misses the path integral measure that is necessary for a sensible 
gravitational theory. Admittedly, the non-perturbative measure of the gravitational path integral is un¬ 
known (although there is numerical evidence that a Lorentz violating preferred foliation might be necessary 
to yield physical results, e.g., [31[). At the perturbative level we know that the role of the measure, along 
with the constraint equations, is to freeze out the non-dynamical degrees of freedom, i.e. scalar and vector 
metric perturbations, in terms of matter degrees of freedom. Therefore, in this admittedly non-covariant 
description, it is not possible to integrate out matter without integrating out part of the metric, which 
consequently leads to either a non-covariant or non-local effective action. In fact, this is exactly what we 
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did in this in paper! Even though well-defined covariant measures exist in simpler gauge theories (such as 
Yang-Mills theory), to our knowledge, one does not exist for quantum gravity. 

From the point of view of canonical quantization, our result shows the non-local structure of the kine¬ 
matic phase space (or Hilbert space) of QFT+gravity, imposed by elliptic GR constraint equations. At an 
intuitive level, this is similar to the Heisenberg uncertainty principle for momentum/position. However, in 
QFT+gravity the uncertainty is between the definition of the vacuum state (or QFT+gravity observables) 
in the UV and the IR. 

In a broader context, the applicability of EFT in quantum gravity is already limited, due to the old 
cosmological constant problem, which suggests extreme fine-tuning of UV physics to explain our observable 
universe on large scales. On a more a theoretical ground, one may consider the firewall paradox [32, 33] 
in the evaporation of black holes as another piece of evidence that local EFT in quantum gravity can fail 
even at very low energies, e.g., at the horizons of large black holes. Indeed, the non-locality introduced by 
imposing the GR constraint equations on the Hilbert space of quantum gravity, as discussed above, has 
been proposed as a possible resolution to the firewall paradox [34]. 
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Appendix A: Kallen-Lehmann Spectral Densities for Stress-Energy 


In this Appendix we describe the tensor structures for stress-energy in the Kallen-Lehmann spectral 
representation, and give the spectral densities /5o,2 in terms of sums over states in Hilbert space. 

Given the general form of Eq. (22), we can extract po or p 2 by contracting respectively with 



(Al) 


(A2) 


Since 



(A3) 


and 


= 0 . 


(A4) 


Consequently, transforming Eq. (22) to momentum space, we find 


Si;‘'‘^^{T^Ak)Tc.p{k'))c,s = {2TT)^S\k + k')2Trpoi-kl), 


(A5) 


and from the second contraction, 


S^''‘^^{T^Ak)Tc.p{k'))a,s = {2TT)^6\k + k')27rp2{-kl). 


(A6) 
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The spectral densities can be expressed in terms of a complete set of states, 

1 = | 0 )( 0 |+ 5 I / + M2)0(po)|p,n)(p,n|. 


(A7) 


Here, p is the total three-momentum for a given state |p,n), and — p^ is its invariant mass- 

squared, with the sum including continuous parameters such as relative momenta of particles. Inserting 
this sum between the operators in {T^i,{x)Tap{ij))s, and transforming to Fourier space as in Eqs. (A5)- 
(A 6 ), we identify the spectral density as 


= k^'^“^( 0 |T^.( 0 )|k,n)(k,n|r„^( 0 )| 0 )+ c.c. 


(A8) 


with p 2 defined similarly with . Only multi-particle states contribute to the sum over intermediate 

states. The contribution from the vacuum cancels when the disconnected part of the correlation function 
{Tfiu){Tai 3 ) is subtracted off, and ( 0 |T^,y|l particle) = 0 because is quadratic in field operators and 
will only act on ( 0 | to create states with two or more quanta. 


Appendix B: An Observational Limit on a Plateau of the CMB Anisotropy Angular Power 

Spectrum 

In Section V C, we showed that the Integrated Sachs-Wolfe (ISW) effect, caused by the metric perturba¬ 
tions that are sourced by vacuum fluctuations add a plateau to the power spectrum of CMB anisotropies 
A 2 ^ ^ Here , we derive an approximate upper limit on this plateau from current CMB observations. 

Given that ACDM predictions (due to Silk damping) as well as measurements of are dropping rapidly 
at high Fs, the strictest bounds on a putative plateau would come from measurements at the highest Vs. 
The best current measurements are provided by the South Pole Telescope (SPT), at the highest Vs (< 
3000) [.jo]. In order to estimate an upper limit on the size of an ISW plateau, we note that missing an 
additional term in the model causes a systematic increase in the x^, given by: 

(xLerved) > x" + [(A?)'®^] ' E ^ ’ (^1) 

I ' 

where the inequality is due to other potential systematic errors (further increasing x^)j while cr; is the 
SPT measurement error for in each 1-bin. x^ iii Eq. (Bl) is a random number which follows the x^ 
distribution: 


p(x^)dx" = 


x” 2 exp (-xV2) rfx^ 


2 "/ 2 r(n/ 2 ) 

n = [# of I — bins] — [# of fitted parameters]. 


(B2) 

(B3) 


Note that we are also ignoring potential degeneracies between the ISW plateau and other cosmological 
parameters that are already fitted for. 

For the ACDM model best fit to SPT data Xobserved = 45.9, which uses 47 Z-bins and 9 fitted parameters. 
Using the errors provided in Table 2 of [-jo], we can plug this into Eq. (Bl), and then integrate the x^ 
distribution to find: 


(B4) 


at 95% confidence level. 


(A2)ISW ^ 2.2 







